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1. INTRODUCTION 

As the lowest "dipole" excitations of the QCD vacuum, 
the light vector mesons (the p meson, in particular) have 
traditionally played an important prototype role in calcu- 
lations and discussions based on QCD sum rules [l[ . In- 
medium versions of these sum rules have been used to set 
constraints on the way in which vector meson masses un- 
dergo possible changes in dense and hot hadronic matter 
0, H, Questions were raised, however, concerning the 
interpretation of such studies. In-medium changes of me- 
son properties, such as their mass shifts in nuclear mat- 
ter, have their primary origin in long-distance physics de- 
scribed by meson-nucleon forward scattering amplitudes 
H and not in the short-distance physics represented by 
subleading terms of the operatorproduct expansion (see 
also related discussions in Refs.jO, In-medium QCD 
sum rules have nonetheless been further developed and 
applied over the years [1, [^, EH, including studies 
with emphasis on the density dependence of four-quark 
condensates [lU, EH ■ The present work aims in a differ- 
ent direction: namely, identifying the spontaneous chiral 
symmetry breaking scale, 4^/^ ~ 1 GcV, and its possi- 
ble change with increasing baryon density, in the context 
of QCD sum rules for the lowest moments of the vector 
meson spectral functions. 

The issue of in-medium changes of hadron properties 
persists as a fundamental theme ever since the Brown- 
Rho (BR) scaling hypothesis [l4| was launched, estab- 
lishing a conceptual relationship between the shifts of 
hadron masses in matter and the sliding scale of spon- 
taneous chiral symmetry breaking with changing ther- 
modynamic conditions. Investigations along these lines 
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included various model calculations of vector meson spec- 
tral functions at finite temperatures and baryon densi- 
ties (see Refs.[l5|, EH, [TtJ and further studies concern- 
ing BR scaling in the context of in-medium QCD sum 
rules, e.g. in Ref . [l8l| ) . Such calculations were performed 
with the aim of understanding the "low-mass enhance- 
ments" observed in dilepton spectra produced in high- 
energy heavy-ion collisions by the CERES/NA45 and 
NA60 [13 experiments at the CERN SPS. These explo- 
rations, primarily focused on the behaviour of the p me- 
son in the strongly interacting hadronic medium, were 
conducted for a long time with two seemingly opposing 
quests: whether there is an in-medium shift of the p me- 
son; or on the other hand, whether the strong collisional 
broadening of the spectral function due to interactions 
of the p meson with nucleons and mesons in the medium 
would render the primary issue of a mass shift physically 
meaningless. 

In the present paper we point out that playing the 
notions of "mass shift" and "broadening" against one 
another may in fact not be the proper question to ask. 
For resonant states such as the p meson, which start out 
with a large decay width already in vacuum, identify- 
ing a mass in an even broader in-medium spectral dis- 
tribution makes sense only in terms of the first moment 
of this spectral distribution. For the two lowest spec- 
tral moments, however, quite accurate statements can 
be made within the framework of QCD sum rules, as we 
shall demonstrate. We propose therefore to abandon the 
"mass shift" versus "broadening" dispute altogether and 
concentrate on an analysis of spectral moments in the 
context of QCD sum rules. Identifying the chiral sym- 
metry breaking scale in such an analysis, both in vac- 
uum and in-medium, permits addressing and examining 
the BR scaling hypothesis in a refined and better focused 
way. 

The strategy pursued in this paper is an update of 
previous work [21| which is in turn closely related to fi- 
nite energy sum rules (FESR) [22|, HH . The advantage of 



these sum rules is that they do not have to rely on the 
existence of a window of stability for the Borcl parameter 
usually employed in the sum rule analysis. Caution must 
nevertheless be exercised with FESR's [24|, [H[ concern- 
ing their sensitivity to high-energy properties of spectral 
functions and the detailed modeling of the transition be- 
tween resonance and continuum regions, a question that 
we shall also address. We concentrate here on the rho 
meson. Starting with vacuum sum rules for the p we re- 
call how the delineation of scales between resonance and 
continuum parts of the spectral function can be related 
to the scale for spontaneous chiral symmetry breaking, 
471-/* ~ 1.2 GeV (the "chiral gap"), where / ff = 92.4 
MeV is the pion decay constant. In- medium sum rules 
are examined using two complementary spectral func- 
tions as generic examples: the one calculated in Ref. [J 
using a chiral meson-nucleon effective Lagrangian with 
vector mesons as explicit degrees of freedom; and the one 
calculated in Ref.[l5j using a model which emphasizes 
the role of particle-hole excitations including baryon res- 
onances. Both types of spectral functions were applied 
earlier fl6l .1281 in descriptions of the CERES/NA45 dilep- 
ton data [19| . Updated versions of such spectral distri- 
butions have been used recently 
with the more accurate NA60 data 
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REMINDER OF QCD SUM RULES FOR 
ISOVECTOR CURRENTS 



We begin with a brief introductory recollection of 
the QCD sum rule approach for excitations carrying 
the quantum numbers of the p meson, J* = 1~ and 
isospin 1 = 1. The corresponding quark current j^(x) = 
^(m7 m u — dj^d) figures in the current-current correlation 
tensor 



U^(q)=i / d 4 x e^ x (Tf(x)f(0)) 



(1) 



In vacuum this tensor can be reduced to a single scalar 
correlation function, H(q 2 ) = \^g llv U. tlv . In a nuclear 
medium the distinction needs to be made between longi- 
tudinal and transverse correlation functions. For vanish- 
ing three-momentum (q^ — (ui,q = 0), the case consid- 
ered here throughout), the longitudinal and transverse 
correlation functions coincide and will again be denoted 
as n(w, q = 0). 

The next step is to write H(q 2 ) as a twice subtracted 
dispersion relation: 

ImII(s) 



IL(q 2 ) = n(0) + II'(O) q 2 + — [ ds- 

7T J 



\ ■ (2) 
s z (s — q — te) 

Alternatively, the same quantity is expressed at large 
spacelike q 2 = — Q 2 < in terms of the Wilson oper- 
ator product expansion (OPE): 



127r 2 IT(g 2 = -Q 2 ) 
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(3) 



In vacuum and for the p meson channel, the expansion 
coefficients are given as: 



Co 

Cl 
C2 



1 



n 



--(m 2 u +m 2 d ) 



(4) 



G 2 ) +6ir 2 (m u (uu) + nid{dd)) 



These three leading coefficients are well determined. The 
dominant perturbative QCD piece Co is shown here in- 
cluding just the standard 0{a s ) correction. At a later 
stage and in all explicit calculations, the QCD correc- 
tions will be further extended up to and including 0{a 3 s ) 
(see Appendix A). 

The quark mass term c\ is small and can safely be 
neglected. The coefficient C2 involves the QCD conden- 
sates of lowest dimension four. The quark condensate 
times the quark mass is given accurately through the 
Gell-Mann - Oakcs - Rcnncr relation as 



(m u uu + rrid dd) ~ m q (uu + dd) 

= ~ m \ f 2 = —(0.11 GeV) 4 



(5) 



The gluon condensate {{u s /n)G 2 ) - (0.3 GeV) 4 is (far 
less accurately) determined by charmonium sum rules. 
For a detailed discussion see Ref. [3~fl | where an upper limit 

((a s /vr) G 2 ) 1 / 4 < 0.31 GeV 

is given. 

In-medium corrections to leading order in the baryon 
density p are introduced by the replacement C2 — > c-i + 
5c 2 {p), with 



5c 2 



3^ 2 [A X M N - —M^ + 2a N ] p 



(6) 



The first term in brackets is the leading density depen- 
dent perturbative QCD correction. It involves the first 
moment, A± = 2(x) u +d, of the parton distribution in the 
nucleon. Given the empirical (MRST) [32|, [33| momen- 
tum fraction carried by u and d quarks in the nucleon, 
(x) u+d ~ 0.62 at Q 2 = 1 GeV 2 , we use A x ~ 1.24 (sec 
Appendix B). 

The second term on the r.h.s. of Eq.© is the correc- 
tion to the gluon condensate at finite density. It is pro- 
portional to the nucleon mass in the chiral limit for which 



= 0.88 GeV from Ref.|34[. The third term 



represents the leading density dependence of the quark 
condensate. It is pro por tional to the nucleon sigma term, 
<jn = (45 ± 8) MeV [35|]. By far the largest contribution 
to #C2 evidently comes from the A\ term, so that the 
large uncertainty in ojss has only relatively minor conse- 
quences. 

Following these considerations the input for C2 and 
<5c 2 is summarized in Table I. The in-medium sum rule 
analysis will be done at normal nuclear matter density, 
p = pa = 0.17 fm~ 3 . 
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value 


reference 


M N 
m q (qq) 
<lfG 2 } 
A 1 
M<™ 
on 


939 MeV 
-(0.11 GeV) 4 
0.005 ± 0.004 GeV 4 
1.237 
0.88 GeV 
45 ± 8 MeV 


GOR 
[31] 
[33] 
[3J 
[35] 



TABLE I: Input summary 



The cocfucicnt C3 involves four-quark condensates in 
the following combination: 

c 3 = -67r 3 a s [((?27 AI 75A a w- J 7A1 7 5 A Q d) 2 ) 

^-{{u lli \ a u + d^X a d) ^ a l)] (?) 

q—u,d,s 
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These condensates of dimension six are not known at any 
reasonable level of precision. What is commonly done at 
this point is to introduce a factorization approximation, 
truncating intermediate states by the QCD ground state 
and writing 



16 
~9 



n(qq)' 



(8) 



with k introduced to parametrize deviations from exact 
factorization (« = 1). The in- medium analogue including 
terms linear in the density p becomes 



448 3 / 2 cr N {qq) \ 

C3 = — —k(p)tt a s [{qq) H p) 

11 \ m„ / 



<?n (qq) 



(9) 



with a density dependent k parameter. 

Clearly, a QCD sum rule analysis that aims for accu- 
racy must try to avoid the uncertain four-quark conden- 
sate piece C3 in the OPE hierarchy. This is indeed possi- 
ble when considering only the two lowest moments of the 
spectral function, ImII(s), as follows. We introduce the 
dimensionless spectral function 



R(s) 



12tt 



ImII(s) . 



(10) 



Note that, in vacuum, R(s) is identified with the ob- 



servable <j(e^ 



hadrons)/cr(e^ 



p + p ). Now 



assume as usual that there exists a delineation scale Sq 
which separates the low-mass resonance region (s < So) 
from the high- mass continuum (s > sq): 



R(s) = R p (s) 9(s -b) + R c (s) 6(s - s ) 



(11) 



This step function delineation between resonance and 
continuum seems schematic on first sight. In prac- 
tice, the transition to the continuum is smooth and sq 
should be considered as an average scale characterizing 
the transition region. A detailed analysis, to be described 
later, shows that the step function ansatz is equivalcntly 



as valid as a more realistic modeling of the threshold 
"ramp", e.g. by the dotted line in Fig[T] 

Let the high-mass continuum be subject to a pertur- 
bative QCD treatment, following duality considerations: 



R c (s) — > c for s > s 



(12) 



Then perform a Borel transformation on Eqs.([2][3]), lead- 
ing to 



127r 2 Il(0) 



ds R(s) e 



= c M 2 + a + 



-s/M 2 



C3 



(13) 



M 2 2M 4 



Choose the (otherwise arbitrary) Borel scale parameter 
sufficiently large, M. > y/so, expand e~ s l M and arrange 
term by term in inverse powers of A4. The result is a 
hierarchy of sum rules for moments of the low-mass part 
of the spectral function R(s): 



/ dsR p (s) = s c + ci - 12tt 2 LT(O) , (14) 
Jo 

/ ds sR p (s) = -£c - c 2 , (15) 
Jo * 

/ dss 2 R p (s) = ^c + c 3 . (16) 
Jo 



These equations are written again to first order in a s , 
with Co = (3/2)(l + a s /ir). Corrections to order a 3 are 
included by the replacements c — > c + (3/2)e n in the 
n-th moment, with e n given explicitly in Appendix A. In 
the detailed calculations the relevant running coupling 
is to be taken as a s (s ) with s ~ 1 GeV 2 , the onset 
scale for the (multipion) continuum part of the quark- 
antiquark excitation spectrum. We use 



a s (s ~ 1 GeV 2 ) = 0.50 ± 0.03 , 



(17) 



referring to the most recent NNLO ( MS ) analysis in 
[36|, The error in ce s (so) is actually the major source 
of uncertainty in the sum rule calculation, all other cor- 
rections being considerably smaller in magnitude relative 
to the leading term. 

The subtraction constant LT(0) in Eq. lfbl)) vanishes in 
vacuum. At finite density this is the Landau term, 
11(0) = jjj— j analogous to the Thomson limit in Comp- 
ton scattering. 

Note that even for a broad spectral distribution R(s), 
a squared "mass" associated with the low-energy sector 
of this spectrum can be defined through the ratio of the 
first and zeroth moments, Eqs. (|14|15p (see also Ref . fiol] ) : 



jo' ds sR(s) 
Io° ds R ( s ) 



(18) 
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3. VACUUM SUM RULES 

A. Identifying the spontaneous chiral symmetry 
breaking scale 

Consider now first the sum rule for the isovector 
current-current correlation function in vacuum. Follow- 
ing Ref. [23[ we start from the working hypothesis that 
the scale sq delineating low-energy and continuum parts 
of the vector-isovector quark-antiquark spectrum should 
be identified with the scale for spontaneous chiral sym- 
metry breaking in QCD: 



LO 



(19) 



For illustration, recall the schematic (large N c ) example 
of a zero-width p meson, 



12tt 2 m 2 n 



(20) 



with the vector coupling constant g. Neglecting small 
quark masses as well as QCD and condensate corrections 
in Eqs. (|14ll5p . one arrives at 



J°daR p (a) = |a = 24^ 2 /' 

/■so 

/ dssR p (s 
Jo 



4 S ° 



(21) 



192^/. , 



and immediately recovers a celebrated current algebra 
result (the KSRF relation HH), 



V2gU 



(22) 



together with the universal vector coupling g = 2-k. 

While this schematic example underlines the validity of 
the hypothesis (fl"9j) . a more detailed test using a realistic 
spectral distribution R{s) and the full sum rule analysis, 
including corrections, must of course be performed. We 
do this along the lines of Ref. [23| and update the results 
found in that work. 

The input is now the resonant p meson spectral func- 
tion R p (s) calculated from one-loop chiral tttt dynamics 
with gauge coupling to vector mesons [3, [39[ . The n-pion 
continuum R c (s) (with n > 4 even) is parametrized as 
in Eq . (fTTj) . with the gap scale sq to be determined by 
the sum rules for the lowest two moments, Eqs. (|14ll5l) . 
The spectral function R(s) is shown in comparison with 
experimental data in Fig[T] 

The analysis proceeds as follows. The equations for 
the two lowest moments of R(s), 



I 



dsRp(s) = s ( c + -£o 



dssR p (s) = 



CO 



2 £l 



+ Ci , 



C2 



(23) 
(24) 



are solved to determine s . For the zeroth moment 
Eq.® gives ^ = 1.13 ± 0.02 GeV. Overall consis- 
tency requires that the same sq results also from Eq. 
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FIG. 1: Vector-isovector spectral function in vacuum showing 
the p resonance and continuum parts as described in the text 
and compared to e + e~ — > tt + tt~ Jo resonance region) and 
e + e~ — > nn data with n even [26l. [27|. 



within an error band determined by the uncertainties of 
the input summarized in table I and Eq. (fT7|) . This test 
turns out to be successful. The detailed analysis of un- 
certainties performed with Eq. (|24p for the first moment 
is shown in FigH The resulting ^ = 1.14 ± 0.01 GeV 
is within 2% of the empirical 4:irf n ~ 1.16 GeV using 
the physical value / w = 92.4 MeV of the pion decay con- 
stant. The postulate (fl9|) identifying ^/so with the scale 
characteristic of spontaneously broken chiral symmetry, 
appears to be working quantitatively. 




1.12 1.14 1.16 1.18 1.20 
^ [GeV] 

FIG. 2: QCD sum rule analysis of the p meson spectral func- 
tion in vacuum. First moment (solid line, left-hand side of 
Eq. (|24|) ') is plotted versus right-hand side (grey band includ- 
ing uncertainties) as function of the gap scale y^so delineating 
low-mass resonance region from high-mass continuum. 

The relation between first and the zeroth moment, 

/•so r s o 

/ dssR p (s) = F(s ) / dsR p {s) (25) 
^o Jo 
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thus involves a uniquely determined function of Sq- 



si (c + - 2c 2 



2s (co 



2ci 



(26) 



up to the estimated uncertainties in the quantities Cj and 
e n (the largest error being associated with a e (s )). The 
squared mass given by fhp = !F(sq) — 0.611±0.013 GeV 2 
or fh p ~ 0.78 ± 0.01 GeV, is very close to the physical p 
meson mass as expected. In fact the canonical relation 
fhp = y/so/2 = s/2 ■ turns out to be satisfied again 
at the 2% level, demonstrating the smallncss of the next- 
to-leading QCD corrections and of the condensate term 

C2- 



B. Sensitivity to continuum threshold modeling 

The question arises whether the quantitatively success- 
ful identification of the continuum threshold yj~s~v with the 
chiral symmetry breaking scale (i.e. the consistency of 
the QCD sum rule analysis with current algebra results) 
is influenced by the schematic step-function parametriza- 
tion (fTTj) . A test can be performed replacing the step 
function by a ramp function to yield a smooth transition 
between resonance and continuum region, as follows: 

R(s) = R p (s) 6(s 2 - s) + R c (s) W(s) , (27) 

where the weight function, W(s), is defined as 



W(x) 



for x < s% 
x — s\ 



for s\ < x < s 2 



S2 - Si 

1 for x > s 2 ■ 



(28) 



The step function behavior is recovered for W(x) in the 
limit Si — > s 2 . 

Using the function W(s), the modified sum rules for 
the lowest two moments of the spectrum R(s) become 



dsR p (s) = s 2 I Co -t- -to 
o \ 2 



ci - 127r 2 n(0) 



-(co-R p (s 2 )) / dsW(s) , (29) 



r ; D ^ 4 ( 3 

J dssR p (s) = — I c + -ei 



<"2 



- (e& - -Rp(sa)) / dssVF(s) . (30) 

J si 

Sets of intervals [si, s 2 ] are then determined so as to sat- 
isfy both sum rules (|29I30[) . and the scale sq defined by 



so 



si + s 2 



respect to variations in the slope (s 2 — si) -1 of the ramp 
function W(s), thus confirming that the step function 
paramctrization of the continuum is not restrictive: the 
smooth "ramping" into the continuum 1 produces values 
of y^So that fall within the narrow (less than 1 %) un- 
certainty band of the step function approach. We note 
at this point that the best fit to the empirical spectral 
function has s 2 — Si — 1 GeV 2 (see FigfT]). It can be 
concluded that the present sum rule analysis and the ob- 
served quantitative agreement of the continuum thresh- 
old with the chiral gap Airf-^ do not depend on details of 
the threshold modeling. 
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FIG. 3: Dependence of ^Tso (determined from Eas. (|29H3T|) ) 
on the slope (s2 — si) -1 of the ramp function W(s) describing 
the onset of the continuum in the vacuum sum rule. The grey 
band indicates the uncertainty range of the result obtained 
with step function parametrization of the continuum. 



4. IN-MEDIUM SUM RULES 

In this section the approach just described is applied 
analogously to vector current spectral functions at finite 
density. We start again from Eqs. (|23l24p . now with in- 
clusion of 11(0) = jjj— and the density dependent correc- 
tions to the condensate terms, c 2 —> c 2 + Sc 2 (see Eq.©). 

Two generic prototypes of in-medium isovector vector 
spectral functions, Imil(w = \fs : q = 0;p), are used for 
demonstration: the one derived from a chiral effective 
Lagrangian with vector meson couplings constrained by 
vector dominance Q (referred to as KKW), and the one 
calculated with emphasis on particle-hole excitations in- 
corporating baryon resonances [l6| (referred to as RW). 
The analysis is performed at the baryon density of nor- 
mal nuclear matter, p = po = 0.17 fm~ 3 . The KKW and 
RW spectral functions, taken at this density, are shown 
in comparison in Figf4j 



(31) 



is now introduced to characterize the continuum thresh- 
old. As shown in Fig[31 the resulting ^/sq is stable with 



1 In this test the uncertainties of a a (Q 2 ) and of the gluon conden- 
sate have been excluded for simplicity. 
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FIG. 4: In-medium isovector vector spectral functions at nu- 
clear matter density, po = 0.16 fm~ 3 , taken from Refs.0] 
(KKW) and HH (RW). The p meson spectrum in vacuum is 
also shown for comparison. 



The KKW and RW in-medium spectral distributions 
both consistently show a strong broadening as compared 
to the vacuum p meson. They differ in details at the 
low mass end of the spectrum. While KKW emphasizes 
the role of chiral in-medium irn interactions. RW focuses 
on the role of nuclcon-holc, A(1232)-hole and TV* (1520)- 
hole excitations. At first sight, none of these broad dis- 
tributions permit identifying an "in-medium mass" or a 
shift thereof with respect to the p meson mass in vacuum. 
This has generally led to the conclusion of there being no 
p mass shift at finite density, but just an overwhelmingly 
large inelastic width due to interactions of the coupled 
p <-> 7T7T system with nucleons in the nuclear medium. 

We now perform the sum rule analysis, first with step 
function continuum, for the two leading moments of the 
KKW and RW spectral distributions: 



ds R p (s) 



s o 1 c o + ^£o 



Cl 



3ir 2 p 
Mn ' 



(32) 



dssR p {s) = ^_ fco + jU ) 



where the gap scale w Sq is permitted to adjust itself to 
the in-medium situation. Consistency of the first and 
zeroth spectral moments is again tested and observed to 
be satisfied within the uncertainties of the input. This 
determines Sq at given density p — po- Effects of smooth 
ramping into the continuum will again be examined later. 

Fig[5] shows the outcome of this procedure for the 
KKW spectral function. In this case, at nuclear mat- 
ter density po, the in- medium gap scale y/s^ is indeed 
seen to be shifted downward from its vacuum position, 
/sb" - 1.14 GeV ~ 471-/^. One finds 



: (1.00 ± 0.02) GeV (KKW at p = po) • (34) 
For comparison, the cross check with the sum rule for the 



zeroth moment gives ^/s^ = (1.02±0.03) GeV, consistent 
with ([341) . 
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FIG. 5: QCD sum rule analysis of the KKW in-medium spec- 
tral function [j. First moment (solid line, left-hand side of 
Eq.(|33}) is plotted versus right-hand side (grey band including 
uncertainties) as function of the in-medium gap scale v/sj. 



(35) 



The analogue of Eg. ([25]) becomes: 

r s o 

ds sR(s, p) = F(sq, p) I dsR(s,p) 



with 



s5 2 (c + §£i) - 2(c 2 + Sc 2 {p)) 
2 [** (c + §e ) + ci - 3ir 2 p/M N ] 



(36) 



The average in-medium "mass" determined from the ra- 
tio ^(sq , p) of the first and zeroth spectral moments is 
found to be 



m*(p) = <JH s hP) = (0.67 ±0.02) GeV (37) 

for the KKW spectral function at density p = po- One 
notes now that the ratio of in-medium and vacuum 1st 
spectral moments behaves as 



m p (vac) 



1 H4i£) 

F(s o ,p = 0) 



0.85 ±0.02 



(38) 



at p = p . 

The successful identification i/s~o = Airf^ in vacuum 
suggests a corresponding generalization to the in-medium 
case: = iirf*, in terms of the pion decay constant, 
/* = ft(p), related to the time component of the axial 
current at finite density. Then one observes \/sq/sq — 
— 0.88 ± 0.02. One finds, within uncertainties, 

_ TO * - ~ & ~ 1- (0.15 ±0.02)-?- , (39) 
m p (vac) Ik po 

suggesting that the BR scaling tendency is indeed visi- 
ble for the KKW in-medium spectral function, contrary 
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to first impression when looking just at the very broad 
overall spectral distribution (3). In this context we refer 
to the subsequent section for an update on the relation- 
ship between the in-medium pion decay constant and the 
density dependence of the chiral condensate. 

The KKW spectrum is based entirely on chiral pion 
dynamics with vector mesons. Baryon resonances are 
assumed to develop large widths and "dissolve" in nu- 
clear matter so that they become part of the continuous 
background. In contrast, the RW spectral function starts 
from a different scenario in which baryon resonances play 
a distinguished role, assuming that they maintain their 
quasiparticlc properties in matter. It is thus instructive 
to conduct, as before, a corresponding sum rule analysis 
for the moments of the RW spectrum under such aspects. 

The result is displayed in FigJFJ] One deduces 



(1.09 ± 0.01) GeV (RW at p = p ) 



(40) 



and y/st/so = 0.97 ± 0.01, together with --f^ ~ 
0.96 ± 0.02 at p = po. (For comparison, the sum rule 
for the zeroth moment gives \/io = (1-11 ± 0.02) GeV, 
consistent with ([40]) ). So the RW spectral function ex- 
hibits dominantly broadening with almost no in-medium 
shift of the ratio of the moments. Notably, both RW 
and KKW based spectral functions work quite well in 
comparison with dilcpton data taken at SPS energies (as- 
suming models for the expansion dynamics of the hot and 
dense matter which have their own uncertainties). This 
implies that it is presumably not possible to distinguish 
between the BR scaling scenario and other (opposing) 
views from those data. 




1.00 
1.04 



1.06 1.08 1.10 1.12 

'4 [GeV] 



1.14 



1.16 



FIG. 6: QCD sum rule analysis of the RW in-medium spec- 
tral function [l6j]. First moment (solid line, left-hand side 
of Eq. (|33|l ) is plotted versus right-hand side (grey band in- 
cluding uncertainties) as function of the in-medium gap scale 



The "ramping" test in order to establish stability with 
respect to the modeling of the continuum is performed 
as for the vacuum case described in the previous section, 
with the same ramping function W(s) employed also for 



the in-medium case. The results of this test for the KKW 
and RW spectral functions are shown in Fig(7] One finds 
again that the determination of y^s^, using a variety of 
smooth transitions to the continuum, is insensitive to 
details of the threshold modeling within the narrow band 
of uncertainties. 



1.03 




Slope of W(s) [GeV ] 

FIG. 7: Dependence of ^/s^, as in Fig|3] on the slope 
(s2 — si)~ of the ramp function W(s), now describing the 
onset of the continuum in the in-medium sum rules. Upper 
panel: result for the KKW spectral function. Lower panel: 
for the RW spectral function. The grey bands indicate the 
uncertainty ranges of the results obtained with step function 
parametrizations of the continuum. 



5. IN-MEDIUM PION DECAY CONSTANT AND 
CHIRAL CONDENSATE: SHORT DIGRESSION 

The present QCD sum rule study asserts that the de- 
lineation between low-energy resonance and high-energy 
continuum parts of the spectral function is related to the 
chiral scale, 47r/ 7r , which acts as an order parameter for 
the spontaneously broken chiral symmetry of the QCD 
vacuum. Its in-medium change with increasing baryon 
density is of fundamental interest and deserves an added 
short section with an updated discussion. 

In the nuclear medium, the relevant quantity is the 
pion decay constant ft(p) = f%(p) related to the time 
component of the axial vector current. Its connection 
with the density dependent chiral (quark) condensate 
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(ipip)p is determined by the in- medium analogue of the 
Gell-Mann, Oakes, Renner relation, 



/:!:_! * 2 



II ■■ 



(41) 



to leading order in the quark mass. Here m* (p) is the 
(charge averaged) pion mass in the medium. A low- 
density theorem gives the leading p dependence of the 
quark condensate as 



(&I>)p = W)o 1 



(42) 



where ctjv = 45 ±8 MeV is the sigma term of the nucleon. 
Assuming that the pion mass is protected by its Gold- 
stone boson nature at low density, we expect to leading 
order in the baryon density: 



i 



p 



i 



p 

6po 



0.83 



(43) 



0.16 fm 3 and taking = 45 MeV for 



at p = po 
orientation. 

A chiral perturbation theory treatment of in-medium 
pion dynamics suggested instead a difference between 
m* and the vacuum pion mass , which translates into 
a stronger density dependence of the pion decay constant, 
ft(p)/U = 1 - (0.26±0.04) / 9/p . On the other hand, the 
charge averaged in-medium pion mass to leading order in 
the baryon density is given by 



™t 2 (p) =m%- T (+) p , 



(44) 



with the isospin-even forward pion-nuclcon amplitude 
T<+) = An(l + m v /M N ) a (+) taken at threshold, 
Empirically [4l|, the corresponding scattering length 



(1.6 ± 1.3) • 10 is compatible with zero. 



This feature derives from a subtle cancellation of non- 
leading terms which cannot be handled accurately in 
baryon chiral perturbation theory Taken as an empir- 
ical constraint, T^^m^) ~ implies m^(p) ~ ra v at 
low density and hence an approximate scaling of /* with 
the square root of the in-medium chiral condensate as 
in Eq. (|43|) . This behavior is actually consistent with the 
observed energy shifts in deeply bound states of pionic 
atoms I42f and related theoretical calculations [43| (see 
also Ref.|44|V 

A recent theoretical study [45[ gives further support to 
these considerations, through a more general derivation 
of ft{p) which does not have to rely on a detailed evalu- 
ation of the charge averaged in-medium pion mass. The 
basic result of Ref. [4a| is 



£(p) = f t (p) = U 



(45) 



where Z and Z* are the wave function renormalization 
factors of the pion in vacuum and in-medium, respec- 
tively. Their ratio is determined by the pion self-energy 
II(w, q, p), as follows: 



J, = i-An M =o,p) 



Using the low-density expression II 
parametrization T^ + '(u) = — tr/v//^ 
rives at 



= -TW p and the 
(jiJ 2 + . . . one ar- 



fZ(p) 

fn 



1 



(47) 



to leading order in the density. With the slope (3 deter- 
mined by the constraint T^ + \lo = m w ) = and assuming 
higher order terms in the expansion of to be small, 
we arrive back at Eq.d43]): f*(p )/U = 0.83±0.03 when 
the admittedly large uncertainty of the nucleon sigma 
term is included. 

Higher order corrections in the density p, calculated 
using in-medium chiral perturbation theory [46| . can be 
expressed in terms of a density dependent effective nu- 
cleon sigma term with a reduced value at normal nuclear 
matter density, a e ^ (po) = (36 ± 9) MeV, leading to a 
3-4% increase of the ratio /*(po)//tt over the value (|4*3")) . 

Notably, the in-medium QCD sum rule analysis assum- 



ing ^/Sq = 47r/* exhibits chiral scaling of this sort for the 
KKW spectral distribution, whereas this is not observed 
for the RW spectral function. 



UJ = 



6. NOTE ON FOUR-QUARK CONDENSATES 

Given spectral functions which consistently satisfy the 
sum rules for the zeroth and first moments, Eqs. (|14ll5|) . 
one can turn to the second moment (|16j) and try to de- 
duce constraints for the four-quark condensate term C3, 
both in vacuum and in-medium. In particular, one can 
discuss deviations from the frequently used factorization 
assumption for those condensates. As mentioned, factor- 
ization means that the intermediate states produced by 
the quark operators entering Eq. ([7]) are truncated by the 
ground state (vacuum) only. Exact factorization means 
K = 1 in Eqs.© and ©. 

When performing the consistency analysis including 
the sum rule (|16[) for the second moment, it turns out 
in all cases that the correction C3 is required to be much 
larger than the value for a factorized four-quark conden- 
sate (with k = 1): factorization proves to be unrealistic 
under any circumstances. For detailed estimates we take 
a value (qq) ~ —(0.2 GeV) 3 and find the following results: 

i) In vacuum, a lower limit k > 4.5 is observed which 
implies strong deviations from factorization. 

ii) For both types of spectral functions (KKW and RW) 
the minimal k required in-medium (typically k > 3) is 
somewhat smaller than in vacuum. 

The range of uncertainty is generally large in all cases, 
with k typically extending from its lower limit up to 
about twice that value. 

One concludes that the four-quark condensates, enter- 
ing the sum rule at the level of the 2nd moment of the 
spectral function, remain basically undetermined. This 
appears to be at variance with reported attempts to con- 
strain such dimension-six condensates from Borel sum 

(46) 

rules for the nucleon [47] . Our findings confirm that the 
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assumption of ground state saturation for four quark con- 
densates should be handled with caution. In the present 
work the sum rules are released from such a dispute by 
restricting procedures to the Oth and 1st moments of the 
spectral distribution for which quantitative statements 
can indeed be made. 



7. SUMMARY AND CONCLUDING REMARKS 

The present work re-emphasizes the usefulness of QCD 
sum rules for moments of spectral functions (or equiva- 
lently, finite energy sum rules), with focus on the p me- 
son both in vacuum and in the nuclear medium. The 
sum rules for the two lowest spectral moments involve 
only the leading (dimension-four) QCD vacuum conden- 
sates as (small) corrections. With inclusion of perturba- 
tive QCD terms up to order a£, these sum rules permit 
an accurate quantitative analysis, unaffected by the large 
uncertainties from condensates of higher dimension (such 
as the four-quark condensates). 

An important scale parameter in this analysis is the 
gap separating low-energy (resonance) and high-energy 
(continuum) regions of the spectral function. For the 
vector-isovector current correlation function, identifying 
this gap with the scale for spontaneous chiral symme- 
try breaking in vacuum, 47r/ 7r , reproduces time-honoured 
current algebra relations and chiral sum rules character- 
istic of low-energy QCD. The corresponding in-medium 
sum rules for the lowest two spectral moments permit to 
address the "mass shift" versus "collisional broadening" 
issue from a new, more quantitative perspective, mean- 
ingful even for broad spectral distributions such as that 
of the p meson at nuclear matter density. Systematic 
tests have been performed to confirm that the conclusions 
drawn from such analysis do not depend on the detailed 
threshold modeling of the transition between resonance 
and continuum parts of the spectral distributions, even 
with strong in-medium broadening. 

Two prototype examples of in-medium rho meson spec- 
tral functions have been examined from this point of view 
in the present paper. Both of these show substantial 
broadening and redistribution of strength into the low- 
mass region, as compared to the vacuum spectrum. The 
sum rule analysis of the lowest spectral moments reveals 
qualitative differences with respect to their Brown-Rho 
(BR) scaling properties. At the same time, both of these 
spectral distributions account quite well for the low-mass 
enhancements observed in dilcpton spectra from high- 
energy nuclear collisions. So one must draw the conclu- 
sion that BR scaling can presumably not be tested in 
such measurements. 

Given the consistency constraints derived from the first 
two sum rules for the spectral moments, one can then 
proceed to the third sum rule equation in this hierarchy 
(involving the second spectral moment and QCD conden- 
sates of dimension six) and discuss limits for the four- 
quark condensates. The outcome of this study demon- 



strates that the frequently used factorization approxima- 
tion for these condensates is questionable under any cir- 
cumstances, both in vacuum and in-medium. 

In summary, we repeat that QCD sum rules for the 
first two moments of vector spectral functions, when com- 
bined with the spontaneous chiral symmetry breaking 
scale of low-energy QCD, permit a quantitatively accu- 
rate analysis in vacuum, consistent with well established 
current algebra relations. The in- medium analogues of 
these sum rules can be used routinely to clarify and clas- 
sify the properties of vector meson spectral functions in 
nuclear matter. An extension to temperature dependent 
sum rules is in progress with special emphasis on the in- 
teresting issue of p — a\ mixing in a thermal pionic heat 
bath. 
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APPENDIX A: QCD CORRECTIONS 

Following Ref. [23[, the expression for the n-th mo- 
ment (with n — 0,1,2) of the spectral distribution in 
the isovector {p meson) channel is written 



s n+1 ( 3 
dss n R p (s) = -9— C + -En 
n+1 \ 2 



(Al) 



+ (-l) B c n+ i - 12tt 2 n(0) (5„ 



The leading perturbative QCD term on the r.h.s. has 



Co 



with 



M -|- Hs. ), The corrections to 0{or s ) are 



,(2) _ 



,(3) - 



1.641 



-10.28 



2.250 
n+1 ' 
11.38 
n + 1 



1.69 



(A2) 



(A3) 



(n+1) 2 
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In applications using (|A1[) the relevant coupling is a s (so) 
with so ~ 1 GeV 2 . In practice we use a s (lGcV 2 ) = 
0.50 ±0.03 fMEzl. 



APPENDIX B: FIRST MOMENT OF QUARK 
DISTRIBUTION 



An accurate value of A\, 



-4, 



2 dxxya 
Jq 



d) , 



(Bl) 



which determines the dominant part of the in-medium 
modifications in our sum rule analysis, is obtained from 
the MRST2001 fits [33j|. In this analysis parton distri- 
butions of the proton are derived from measurements 
of structure functions by the HI and ZEUS collabora- 
tions at HERA, and by the DO and CDF collabora- 
tions at the Tevatron, performing DGLAP evolution. 
The parametrization of the parton distributions at Q 2 = 



1 GeV is: 



xu v = 


0.158 x 25 (l - 


xd v = 


0.040 x°- 27 (l - 


xS = 


0.222x-° 26 (l - 


xA = 


x(d — u) 




1.195 x 1;24 (l- 


2u 


0.45- A , 


lit = 


0.45 + A , 



)3 .33 ( 
\3.88/ 



""(1 



x) 7 - 10 (l 



5.61a; - 5 + 55.49a;) , 
52.73a; - 5 + 30.65a;) , 
+ 3.42a; - 5 + 10.30a;) , 



\9.10 



(1 + 14.05a; - 45.52a; 2 ) , 



(B2) 



where u v and d v denote the valence u- and d-quark distri- 
butions while 2u and 2d are the sea quark distributions. 
A denotes the difference between d and u. 

Using this parametrization, Ai at a 1 GeV scale is 
directly calculated as 

yli = 2/ dxx(u v +d v + 2u + 2d) = 1.2373 . (B3) 
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